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Abstract We study incremental wave propagation for what is seemingly the sim- 
plest boundary value problem, namely that constitued by the plane interface of a 
semi-infinite solid. With a view to model loaded elastomers and soft tissues, we 
focus on incompressible solids, subjected to large homogeneous static deformations. 
The resulting strain-induced anisotropy complicates matters for the incremental 
boundary value problem, but we transpose and take advantage of powerful tech- 
niques and results from the linear anisotropic elastodynamics theory. In particular 
we cover several situations where fully explicit secular equations can be derived, 
including Rayleigh and Stoneley waves in principal directions, and Rayleigh waves 
polarized in a principal plane or propagating in any direction in a principal plane. 
We also discuss the merits of polynomial secular equations with respect to more 
robust, but less transparent, exact secular equations. 

1 Introduction 

The term "acousto-elastic effect" describes the interplay between the static deformation 
of an elastic solid and the motion of an clastic wave. If both the deformation and the 
motion are of infinitesimal amplitude, then all the governing equations are linearized, see 
for instance the Chapter by Norris for examples and applications or the experimental 
results of Pao et al. (1984). If both the deformation and the motion are of finite ampli- 
tude, then the resulting governing equations are highly nonlinear, and their resolution is 
the subject of much research, see the Chapter by Fu for the weakly nonlinear theory and 
the Chapter by Saccomandi for the fully nonlinear theory. 

In between those two situations lies the theory of "small-on-large" , also known as the 
theory of "incremental" motions, where the wave is an infinitesimal perturbation super- 
imposed onto the large static homogeneous deformation of a generic hyperelastic solid. 
There, the homogeneous character of the static deformation and the linear character of 
the incremental equations of motion ensure that the calculations are valid for any strain 
energy density (to be specified later for applications, if necessary). The next Section 
of this Chapter briefly recalls the governing equations of incremental motions (see the 
Chapter by Ogden for their derivation). 

It turns out that many similarities can be drawn between the equations of incremental 
motions and those of linear anisotropic elasticity, with the main difference that in the 
latter case, the anisotropy is set once and for all for a given crystal whereas in the former 
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case, it is strain-induced and susceptible to great variations from one configuration to 
another. Using the similarities, we may transpose the so-called Stroh formulation and 
exploit its many results; on the other hand, when focussing on the differences, we may 
highlight the influences of the pre-stress and of the choice of a strain-energy density 
on the propagation of waves. In this Chapter, attention is restricted to waves at the 
interface of pre-deformed, semi-infinite solids, in contact either with vacuum (Raylcigh 
waves) or with another solid (Stoneley waves). With a view to model elastomers and 
biological soft tissues, the solids are considered to be incompressible (mathematically, 
this internal constraint lightens somewhat the expressions but does not prove essential 
to the resolution). 

Several situations are treated: principal wave propagation in Section 3, principal po- 
larization in Section 4, and principal plane propagation in Section 5. The emphasis is 
on deriving explicit secular equations in polynomial form, using some simple "funda- 
mental equations" derived at the end of Section 2. Of course, as the setting gets more 
and more involved, so does the search for a polynomial secular equation; eventually its 
degree becomes too high for comfort and other techniques are required. The concluding 
section (Section 6) discusses the pros and cons of such equations, as opposed to exact, 
non-explicit, secular equations, free of spurious roots. 

2 Basic equations 

2.1 Finite deformation 

Consider an isotropic, incompressible, hyperelastic solid at rest, characterized by a 
mass density p and a strain energy function W. Then subject it to a large, static, ho- 
mogeneous deformation ( "the pre-strain" ) carrying the particle at X in the undeformed 
configuration to the position x in the deformed configuration. 

Call F — dx/dX the corresponding constant deformation gradient and B = FF* the 
associated left Cauchy-Green strain tensor. This tensor being symmetric, the directions 
of its eigenvectors are orthogonal; they are called the principal axes of pre-strain or in 
short, the principal axes. Also, the eigenvalues of B are positive, Af, A|, A3, say, and 
Ai, A2, A3 are called the principal stretches. Figure 1 shows how a unit cube with edges 
aligned with the principal axes, is transformed by the pre-strain. Note that because the 
solid is incompressible, its volume is preserved through any deformation so that here, 

AiA 2 A 3 = 1. (2.1) 

The first two principal invariants of strain arc defined as 

h = tr B, I 2 = [(tr B) 2 - tr (B 2 )]/2. (2.2) 

In the Cartesian coordinate system aligned with the principal axes, B is diagonal. Calling 
ei, e2, e3, the unit vectors in the x\, X2, £3 directions, respectively, we have 

B = Aid <g> ei + \%e 2 <g> e 2 + A^ <g> e 3 , (2.3) 

and the computation of I\, I2 there gives 

h = \\ + \l + \l h = \l\l + \l\l + \l\l (2.4) 
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Figure 1. Finite homogeneous deformation of a unit cube. 



For an isotropic solid, W may be given as a function of the invariants: W — W(Ii, 1 2) 
or cquivalently, as a symmetric function of the principal stretches: W = W(X±, A2, A3), 
according to what is most convenient for the analysis or according to how W has been 
determined experimentally. With the first choice, the constant Cauchy stress ( "the pre- 
stress" ) necessary to maintain the solid in its state of finite homogeneous deformation 
is 

tr = -pl+ 2{dW/dh + hdW/dh)B - 2(dW/dI 2 )B 2 , (2.5) 

where p is a Lagrange multiplier due to the constraint of incompressibility (a yet arbitrary 
constant scalar to be determined from initial and boundary conditions.) With the second 
choice, the non-zero components of cr relative to the principal axes are written as 

cr l = -p + XidW/dXu i = 1, 2, 3 (no sum). (2.6) 

The proof for the equivalence between (2.5) and (2.6) relies on the connections (2.4). 

2.2 Incremental equations 

Consider a half-space filled with an incompressible hyperelastic solid subject to a large 
homogeneous deformation. We take the Cartesian coordinate system (x\, X2, X3) to be 
oriented so that the boundary is at £2 — 0, and we study the propagation in the x\ 
direction of an infinitesimal interface wave in the solid. 

This wave is inhomogeneous as it progresses in an harmonic manner in a direction 
lying in the interface, while its amplitude decays with distance from the boundary. 

We call u the mechanical displacement associated with the wave, and p the increment 
in the Lagrange multiplier p due to incompressibility. The incremental nominal stress 
tensor s has components 

= AojiikUk,i +puj,i - pStj, (2.7) 
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where the comma denotes partial differentiation with respect to the coordinates Xi, x 2 , 
£3. Here, „4o is the fourth-order tensor of instantaneous elastic moduli, with components 

d 2 W 

A juk = Fj a Fip— — — = Aoikji- (2-8) 

-Fia-Fkp 

Note that due to the symmetry above, Ao has in general 45 independent components. 
In the principal axes coordinate system (xiX 2 Xs) however, there are only 15 independent 
non-zero components; they are (Ogden, 2001): 

Aotijj — A^ Xj ^^ij , 

Aoijij = (XiWi - X j W j )X 2 i /(X 2 i - A, 2 ), i ? j, A, ? X 3 , 

•Aoijij = i-Aoim — Aoujj + AjWj)/2, i ^ j, A« = Xj, 

■Aoijji = Aojuj = Aoijij — AjWj, i ^ j, (2-9) 

(no sums on repeated indexes here), where Wj — dW/dXj and Wij = d 2 W / {dXidXj) . 

Finally, the governing equations are the incremental equations of motion and the 
incremental constraint of incompressibility; they read 

Sjt.j = pd 2 Ui/dt 2 , Ujj = 0, (2-10) 

respectively. 

Now everything is in place to solve an interface wave problem. We take u and p in 
the form 

{u,p} = {U{kx 2 ),ikP{kx 2 )}c ik< - n -^ vt \ (2.11) 

where k is the wave number, U and P are functions of the variable kx 2 only, n is the 
unit vector in the direction of propagation, and v is the speed. Clearly by (2.7), s has a 
similar form, say 

s = ikS(kx 2 )e ik{n - & ' vt) , (2.12) 

where S is a function of kx 2 only. 

After substitution, it turns out that the incremental governing equations (2.10) can 
be cast as the following first-order differential system (Chadwick, 1997), 

£' = \N£, where i = [U,t] t , (2.13) 

the prime denotes differentiation with respect to the variable kx 2 , and t are the tractions 
acting on planes parallel to the boundary, with components tj = S 2 j. Here the matrix 
N has the following block structure 



N = 



JVi N 2 
N 3 + pv 2 I N\ 



(2.14) 



where TVi, AT 2 = N 2 , and N3 = N 3 are square matrices. This is the so-called Str oh for- 
mulation. In effect, many of the results established thanks to the Stroh (1962) formalism 
in linear anisotropic elasticity can formally be carried over to the context of incremental 
dynamics in nonlinear elasticity, as shown by Chadwick and Jarvis (1979a), Chadwick 
(1997), and Fu (2005a,b). 
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2.3 Resolution 

The solution to the first-order differential system (2.13) is an exponential function in 

kx 2 , 

£(kx 2 ) = e ik ^C, (2-15) 

where C is a constant vector and q is a scalar. Then the following eigenvalue problem 
emerges: AT£ = q£. Its resolution is in two steps. 

First, find the eigenvalues by solving the propagation condition, 

det (N - ql) = 0, (2.16) 

for q, and keep those qj's which satisfy the decay condition. For instance, when the solid 
fills up the x 2 ^ half-space, the decay condition is 

> 0, (2.17) 

ensuring that the solution (2.15) is localized near the interface and vanishes away from 
it. The penetration depth of the interface wave is clearly related to the magnitude of 
3(g): the smaller this quantity is, the deeper the wave penetrates into the solid. 

The propagation condition (2.16) is a polynomial in q with real coefficients and it 
has only complex roots (Fu, 2005b), which come therefore in pairs of complex conjugate 
quantities. Hence, half of all the roots to the propagation condition qualify as satisfying 
the decay condition. Let C be the eigenvector corresponding to the qualifying root qj. 

Now proceed to the second step, which is to construct the general localized solution 
to the equations of motion, as 

i(kx 2 ) = £ 7j e ifc «^ 2 C J , (2-18) 

for some arbitrary constants jj. Then compute this vector at the interface x 2 = and 
apply the boundary conditions. The vector £(0) is often decomposed as follows, 



U(0) 
t(0) 



3 — 



7, (2.19) 



where A and B are square matrices and 7 is the vector with components 7^. For instance, 
the archetype of interface waves is the Rayleigh (1885) surface wave, which propagates 
at the interface between a solid half-space and the vacuum, leaving the boundary free of 
tractions. Mathematically, the corresponding boundary condition is that t(0) = -B7 = 0, 
leading to 

det B = 0. (2.20) 

This (complex) form of the secular equation is however not the optimal form, and it might 
lead to unsatisfactory answers to the questions of existence and uniqueness of the wave 
(see Barnett (2000) for an historical account of this point). From the Stroh formalism, 
and its application to the present context, we learn that it is much more efficient to 
work with the surface impedance matrix (Ingebrigsten and Tonning, 1969) than with the 
matrices A and B; this matrix M is defined by 

M = \BA 1 . (2.21) 
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It is Hermitian (Barnett and Lothe, 1985; Pu, 2005b) and so det M = -i(det B)/(det A) 
is a real quantity and the secular equation for Rayleigh surface waves, written in the form 

det M = 0, (2.22) 

is a real equation, in contrast to (2.20). Moreover, if there is a root to this equation in 
the subsonic regime (where v is less than the speed of any bulk wave), then it is unique; 
also, the existence of a root is equivalent to the existence of a surface wave. 

Similar results also exist for other types of interface waves as seen in the course of this 
Chapter. For instance, the boundary conditions for Stoneley (1924) interface waves are 
that displacements and tractions are continuous across the boundary between two rigidly 
bonded semi-infinite solids. Then £(0) = £*(0) where the asterisk refers to quantities 
for the solid in 12 0. Equivalently, A-y — A*j*, B7 = B*j*, from which comes 
[BA~ l A* - B*}j* = 0, leading to 

det(M + M*) = 0, (2.23) 

the optimal form of the secular equation for Stoneley interface waves. Here M* is the 
surface impedance matrix for the solid in the £2 ^ half-space, defined as 

M* =iB*{A*)~ 1 . (2.24) 

2.4 Explicit secular equations 

The derivation of a secular equation, preferably in the optimal form involving the 
surface impedance matrix, is no sinecure in general. The problematic step lies in the 
resolution of the propagation condition (2.16). 

For principal wave propagation (xi, x 2 , £3 are aligned with the principal axes), the 
propagation condition factorizes into the product of a term linear in q 2 and a term 
quadratic in q 2 . Here we can compute the roots explicitly, keep the qualifying ones (see 
(2.17)), and solve the boundary value problem in its entirety. Many problems falling in 
this category have been solved over the years, and some are presented in Section 3. 

For a non-principal wave with propagation direction and attenuation direction both 
in a principal plane (the saggital plane {xix?) is a principal plane but x\ is not a principal 
axis), the propagation condition factorizes into the product of a term linear in q 2 and a 
term quartic in q. We treat this case in Section 4. Although it is possible to write down 
formally the qualifying roots of the quartic (Fu, 2005a; Destrade and Fu, 2006; Fu and 
Brookes, 2006), the formulas involved are cumbersome to interpret. 

For a wave propagating in a principal plane but not in a principal direction (£2 is 
aligned with a principal axis but neither x\ nor x 3 are aligned with principal axes), the 
propagation condition is a cubic in q 2 . We treat this case in Section 5. Now it is a 
daunting task to find analytical expressions for the roots q satisfying the decay condition 
(2.17). 

Finally, for wave propagation in any other case, the propagation condition is a sextic 
in q, unsolvable analytically according to Galois theory. 

These observations suggest that, except in the case of principal waves, numerical 
procedures are required in order to make progress. It is indeed the case that sophisticated 
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tools and efficient numerical recipes have been developed by Barnett and Lothe (1985), 
Fu and Mielke (2002), and several others, with most satisfying results. However it is 
also the case that some interface wave problems can be solved analytically, up to the 
derivation of the secular equation in explicit polynomial form. The first steps in that 
direction were taken by Currie (1979), and his advances were later refined by Taylor and 
Currie (1981) and Taziev (1989), revisited by Mozhaev (1995) and by Ting (2004), and 
extended by Destrade (2003). 

The equations that turn out to be fundamental in the derivation of explicit polynomial 
secular equations are 

£(0) • IN n £(0) = 0, (2.25) 



where 7 : 



7 
7 



and n is an integer. Their derivation is most simple. First, it can be 



shown by induction (Ting, 2004) that N n has a block structure similar to that of N, 
that is 



TV™ = 



N. 

K (n) jyC 



{n) 



(n)t 



with K (n) = K (n) \ N [ 2 n) = N 2 n)t . It then follows that 



IN = 



K {n) jyO 

n: 



(n)i 



(n) 



(2.26) 



(2.27) 



is symmetric for all n. Now take the scalar product of both sides of the governing 
equation (2.13) by 7iV™£ to get 

£' • IN n £ = i£ • IN n+1 £; (2.28) 

finally add its complex conjugate to this equality to end up with 

£' • 77V"f + £ • IN n t = 0, (2.29) 

and, by integration between the interface (at £2 — 0) and infinity (where U and t, and 
thus £, vanish), arrive at (2.25). 

For instance, the boundary condition for Rayleigh surface waves is that there are 
no incremental tractions at the interface; thus £(0) = [17(0), 0]*, and the fundamental 
equations (2.25) reduce to 

U(0)-K {n) U(0) = 0. (2.30) 



3 Principal waves 

Here we take (xi,x 2 ,x 3 ) to coincide with the principal axes (x\, X2, #3). The pre- 
deformation is thus 

xi = X1X1, x 2 = X2X2, x 3 = X 3 X 3 . (3.1) 

Figure 2 summarizes the situation with respect to the waves' characteristics near the 
interface. Bear in mind that the wave analysis is linear and gives no indication about 
the amplitude; moreover, a half-space has no characteristic length so that the secular 
equation is non-dispersive and the wavelength remains undetermined. 
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3.1 Governing equations 

For principal waves, the fields (2.11) and (2.12) are independent of x 3 = x 3 , because 
x 2 — X2 and n • x = X\ = x\. Also, recall from (2.9) that the non-zero components of 
Aq in the {x\,x 2 ,x 3 ) coordinate system of principal axes are 



.4.01111 — AjWil, -4oil22 — A1A2W12, Awn.12 — \ 2 W 22 

r 2- _,-2 A _ XiW 1 -X 2 W 2 _ x 2 w 1 -\ 1 w 2 ^ x 

\ -^01212 — ^2 ^02121 — T2 T2 7 -^01221 — 72 T2 M-*>2, (>J- Z ) 

A^ — \ 2 A^ A 



and also .403333, -4oii33i ^02233^ -4oi3i3j ^023231 -4o3i3i7 ^03232, ^4oi33i, an d ^02332, 
whose expressions are not needed in this Section. 

From these observations follows that the third equation of motion (2.10)3 reduces to 

-S 13 +iS' 23 = -pv 2 U 3 , (3.3) 

where by (2.7), 

LS13 = 14.01313^3, i#23 = ^02323^3- ( 3 - 4 ) 

Hence the movement along the x 3 principal axis is governed by an equation which depends 
only on U3. For this equation, governing what is termed the anti-plane motion, we 
take the trivial solution: U3 — 0, and we focus on the in-plane motion. According to 
(2.10)1,2,4, it is governed by 

-S n +iS' 21 = -pv 2 U 1: -S 12 + iS 22 = - P v 2 U 2 , tf7i + £/^ = 0, (3.5) 
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where by (2.7), 



iS n = i(Aim +p)Ui + A mi22 U' 2 - P, 

-^02121 U[ +i(A)1221 +P)U 2 , 
iSu = (Ami +P)U[ + L4oi212^2, 
1522 = iAl!22^1 + (-402222 +P)U 2 - P. 



(3.6) 



We eliminate p in favour of the pre-stress: by (2.6) at j = 2, we have p = X 2 W 2 — a 2 and 
so by (3.2), 

A1221 +P = A2121 - 0-2- (3.7) 



It then follows from the second equation above that 
U[=i 



-4o2121 — C 2 1 

1^2 H Z iJ21 



-4(32121 



-4(32121 



(3.8) 



and this constitutes the first line of the first-order system (2.13). The second line comes 
from the incremental incompressibility constraint (3.5)3 as 



U^ = i[-U!}. 



(3.9) 



Proceeding similarly for S' 21 , S' 22 , we find eventually that the governing equations are 
indeed in the form (2.13), where £ = [Ui, U 2 , S 2 i, £22]* an d —Ni, N 2 , and — N 3 are 
given by 



q 721 — °2 




r 1 

— 




2(/3i2 +721 - cr 2 ) 


721 




721 







1 










712 





(721 ~ 02) 2 
721 



(3.10) 



respectively, where we used the following short-hand notations (no sums), 



20a = •A-OUii + -Aojjjj 2A.QUjj 2A.Qijji — 2(3 ji 



(3.11) 



or cquivalently, 



Hi = (AiWi - XjW^XViXj - X 2 ) = x 2 xj 2 ljl7 

2^ = X\Wa - 2X i X j W ij + X^Wjj + 2(X i W j - X^X^/iX 2 - X*) = 2f3 M . (3.12) 



3.2 Resolution 

The propagation condition (2.16) reduces to a quadratic in q 2 , 

72i9 4 + (2/3i2 - pv 2 )q 2 + 712 - pv 2 = 0. 



(3.13) 



Notice how a 2 , though present in AT, does not appear explicitly in this equation. 
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Calling q\, q 2 , the roots of the quadratic, we have 



2 2 
Qll2 



712 - pv 

721 



2 1 2 
Qi + ?2 



721 



(3.14) 



The roots qi, q 2 of the biquadratic satisfying the decay condition (2.17) are in one of 
the two following forms; either: q\ = i/3i, 92 = i/?2, where /?i > 0, /?2 > 0, or: q\ = a + i/3, 
(J2 = —a + where /3 > 0. Whatever the case, q\q\ > 0, q\q 2 < 0, and qi + q 2 is 
a purely imaginary quantity. From the first inequality we deduce that (Dowaikh and 
Ogden, 1990) 



V 



1 712 - ptr 
721 



(3.15) 



is a real quantity. From the second, and using the definitions of 77 and 7^ , we find 

./^12 r... 2 



/ , \2 712 - 2/3l2 2 

9i92 = -r?, (<?i+92) = 277 - r? 

721 



A X A 2 — 2 



721 



2??-r/ 2 . (3.16) 



We compute the eigenvectors C 1 and £ 2 of TV corresponding to qi and q 2 as any 
column of the matrix adjoint to TV — q\I and to TV — 52-?, respectively. Choosing the 
third column, we find 



a 1 ' 




a 2 ' 


b\ 


, c 2 = 


b 2 



where 



721 721 



21* 



6 J = [-q 3 (q 2 - 1 + o 2 ), <? 2 (1 - o 2 ) - ry 2 ] 



(3.17) 



(3.18) 



and cf 2 = 02/721 is a non-dimensional measure of the pre-stress. 

We can now construct the A = [a^a 2 ] and B = [b 1 ^ 2 ] matrices, and the surface 
impedance matrix M = — \BA 1 . It turns out to be 



M = — i72i 



Ql + 92 1-CT2-7? 

-(l-o 2 -7?) (qi+q2)v 



(3.19) 



which is indeed Hermitian because q\ + q 2 is a purely imaginary quantity and 77 is real. 
For Rayleigh surface waves, the secular equation is (2.22), or here, using (3.16), 



V + (2 - A 2 A 2 - 2 + 2^ - 2o 2 )77 - (1 - o 2 ) 2 = 0. 
721 



(3.20) 



Dowaikh and Ogden (1990) established this form of the secular equation for principal 
surface waves in pre-stressed incompressible solids, following other works by Hayes and 
Rivlin (1961), Flavin (1963), Willson (1973a,b), Chadwick and Jarvis (1979a), Guz (2002, 
review with an extensive bibliography), and many others. It is of course consistent with 
Lord Rayleigh's own analysis of surface waves in linear isotropic incompressible solids. 
To check this, let the solid be un-stressed (oj = 0) and un-deformed (A, = 1); then r] 
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reduces to yl — pv 2 / /Jo, where po is the infinitesimal shear modulus; also, /3\ 2 — 721 
and 77 is the real root of rf + rj 2 + 3i] — 1 = 0, that is 77 ~ 0.2956 giving pv 2 /p — 0.9126, 
as found by Rayleigh (1885). 

For Stoneley interface waves, the secular equation is (2.23) where 



M + M* = -i 



721(91 + 92) - 72l(9* + Q2) 72i(l - 77) - 7^(1 - T]*) 
-72i(l - J?) +7|i(l - V*) 72i(9i + Q2)v - 7 2 *i ill + <?2) 7 7* 



(3.21) 



This equation was studied in great detail by Dowaikh and Ogdcn (1991) and by Chadwick 
(1995). It is consistent with the analysis of Stoneley (1924) of interface waves in linear 
isotropic incompressible solids. A remarkable feature of this secular equation for principal 
Stoneley interface waves in deformed incompressible solids - first noted by Chadwick and 
Jarvis (1979b) - is that the pre-stress a 2 does not appear explicitly in it, in contrast to 
the equation for surface waves (3.20). This quantity, which is continuous across the 
interface (a 2 = a 2 ), disappears in the addition of the two surface impedance matrices. 
Of course it still plays an implicit role, in determining the pre-strain. 

Dowaikh and Ogdcn (1990, 1991), Chadwick (1995), and Guz (2002, review) have 
covered almost every aspect of principal interface wave propagation and more information 
can be found in their respective articles. In the next Subsection we rapidly work out two 
examples of surface waves. 



3.3 Examples 

First we present an example taken from the literature on elastomers, where the 
Mooney-Rivlin strain energy function is often encountered. It is given by 

W = ©i(A? + A 2 + A 2 - 3)/2 + 2? 2 (A 2 A 2 + A 2 , A3 + A 2 , A 2 - 3)/2, (3.22) 

where T>\ and T> 2 are positive constants with the dimensions of a stiffness. The Mooney- 
Rivlin material enjoys special properties with respect to wave propagation (the neo- 
Hookcan material, which corresponds to the special case D 2 = 0, enjoys even more 
special properties as is seen in Section 5.3). For instance, once subjected to a large 
homogeneous pre-strain, it permits the propagation of bulk waves in every direction; 
these waves can be infinitesimal, but also of arbitrary finite amplitude (Boulanger and 
Hayes, 1992); they can be homogeneous plane waves but also inhomogeneous plane waves 
(Destrade, 2000, 2002). The quantities (3.12) are also quite special; they are 

7ij = (V, + V 2 \ 2 )\ 2 , 2/3y - (©! + V 2 X 2 )(X 2 + \% (3.23) 

where k 7^ and thus they satisfy 

23 i:j = 7,-,- + 7,, : . (3.24) 

These relationships mean that the biquadratic (3.13) factorizes to 

(q 2 + l)(q 2 + n 2 )=Q, (3.25) 

and that the secular equation (3.20) reduces to 

77 3 + 77 2 + (3 - 2a 2 )77-(l-a 2 ) 2 = 0. (3.26) 
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Hence one qualifying root is q\ = i, whatever the values of the material constants T>\ 
and T>2- The other root is q 2 = ir/. When there is no pre-stress normal to the boundary 
(a 2 = 0), then r\ is the real root of rf + rj 2 + 3r? — 1 = 0, that is 77 ~ 0.2956 giving 

pv 2 = 7l2 - 721?? 2 = {T> x + V 2 X 2 )(X 2 - 0.0874A 2 ), (3.27) 

a result first established by Flavin (1963). Here q\ — i and q 2 ~ 0.295i so that the 
penetration depth of the surface wave is fixed and is completely independent of the pre- 
strain and of the material parameters T>\ and T> 2 . We say that the penetration depth is 
universal relative to the class of Mooney-Rivlin materials. 

The second example is taken from the biomechanics literature. From a series of 
uniaxial tests on human aortic aneurysms, Raghavan and Vorp (2000) deduced that the 
following strain energy density gave a satisfying fit with the data plots, 

W = Ci(A? + A 2 + A 2 - 3) + C 2 (A 2 + A 2 + A 2 - 3) 2 , (3.28) 

where, typically, C\ = 0.175 MPa, C 2 = 1.9 MPa (Karduna et al. (1997) use the same 
expression to model the response of passive myocardium.) A uniaxial pre-stress is u\ 7^ 0, 
C2 = C3 = 0, leading through (2.6) to the following equi-biaxial pre-strain, 

Ai = A, A 2 = \- 1/2 , A 3 = A~ 1/2 , (3.29) 

where A is calculated from 

(7i = 2(A 2 - A _1 )[Ci + 2C 2 (2 - 3A- 1 + A~ 3 )]. (3.30) 

Finally using the following expressions for the relevant moduli, 

721 = 2C 1 A~ 1 + 4C 2 (2A- X - 3A~ 2 + A~ 4 ), 7l2 = A 3 721 , 

P21 = Ci(A 2 + A" 1 ) + 2C 2 (4A 2 - 3A - A" 1 - 3A" 2 + 3A~ 4 ), (3.31) 

it is a simple matter to solve the secular equation (3.20) numerically and plot the vari- 
ations of the squared wave speed, scaled with respect to the squared bulk wave speed 
712 / P, with the pre-stretch A. Figure 36 displays these variations; for comparison pur- 
poses, Figure 3a shows the variations of the scaled squared wave speed in the case of a 
Mooney-Rivlin material in uniaxial stress; in that later case the graph is independent of 
the material parameters T>\ and T> 2 because by (3.27), pv 2 /-f 2 i = 1 — 0.0874A~ 3 . The 
dashed lines indicate the speed of Lord Rayleigh's squared speed in the isotropic (no 
pre-strain) case where A = 1, pv 2 /712 — 0.9126. Notice how different the responses of 
two solids are in that neighbourhood. Note also that for high compressive stretches, the 
squared speeds eventually falls off to zero; this happens at A ~ 0.0874 1 / 3 ~ 0.444 for all 
Mooney-Rivlin materials, as shown by Biot (1963), and at A ~ 0.315 for the soft biolog- 
ical tissue model above. Beyond that critical compression stretch, v 2 < 0, leading to a 
purely imaginary v, an amplitude which then grows exponentially with time according 
to (2.11), and a breakdown of the linearized analysis. The search for critical compression 
stretches is an extremely active area of research, clearly linked to the geometric stability 
analysis of solids. 
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Figure 3. Variations of the scaled squared surface wave speed with the stretch in uniaxial 
pre-stress (a) for any Mooney-Rivlin material and (b) for the solid with strain energy 
density (3.28) where C x = 0.175 MPa, C 2 = 1.9 MPa. 



4 Waves polarized in a principal plane 

In this Section we study the case where £3 is aligned with the principal axis X3 but neither 
X\ (propagation direction) nor x 2 (attenuation direction) are aligned with principal axes. 
The components Aojuk in the (X1X2X3) coordinate system of the instantaneous elastic 
moduli tensor Ao are related to the components Aojuk, given by (3.2), in the (X1X2X3) 
coordinate system of principal axes through the tensor transformations 



A 



Qjilk 



pQiq^lrQ ks A()pqr s j 



where Jl 



cos 9 — sin 6 
sin 9 cos 9 
1 



(4.1) 



and 9 is the angle between xi and x\. 

In particular we find that the non-zero components in the forms -Aoi3Zfc and *4o23/fc 
are A1313, A1323, A1331, A2323, A2313, and A 2 332- The mechanical fields (2.11) and 
(2.12) are independent of £3 = X3 because n • x — X\ here. As a consequence, the third 
equation of motion (2.10)3 reduces to 



-S 



13 



■iS'o 



23 



-pv 2 U 3 , 



where by (2.7), 



13 



LA 



01313C3 



U3 + ^01323^3- iSi 



'23 



02313^3 



-402323^3- 



(4.2) 



(4.3) 



Hence the movement along the X3 principal axis is governed by an equation which depends 
only on U3, and for this anti-plane motion we take the trivial solution: U3 = 0. 

The equations governing the in-plane motion have been derived in the case of a 
general plane pre-strain by Fu (2005a) and solved for surface waves in the case of a pre- 
strain consisting in a triaxial stretch followed by a simple shear by Destrade and 9gden 
(2005). Instead of treating these cases again, we revisit the case relative to one of the 
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most important pre-strain fitting into the present context, that of finite simple shear, 
presented originally by Connor and Ogden (1995) for surface waves. 

Figure 4 sketches what happens to a unit cube when a solid is subject to the simple 
shear of amount K , 

x 1 =X 1 +KX 2 , x 2 = X 2l x 3 =X 3 . (4.4) 

Here the principal axes are x 3 = X 3 and x±, x 2 which make an angle if) with X\ and 
with X 2l respectively. That angle, and the corresponding principal stretches are (e.g. 
Chadwick (1976)), 

= (l/2)tan- 1 (2//sT), Ai, 2 = y/l + K 2 /A ± K/2, A 3 = 1. (4.5) 

These relations highlight a major difference between this homogeneous pre-strain and 
the triaxial pre-stretch (3.1): here the orientation of the principal axes with respect to 
the plane interface changes as the magnitude of the pre-strain changes. 




Figure 4. Finite simple shear of amount K of a block near the interface. The x\ direction 
is the direction of shear; the (xi£ 2 ) plane is the plane of shear; the (xix 3 ) plane is the 
glide plane. 



4.1 Governing equations for simple shear pre-strain 

The deformation (4.4) is an example of plane strain. As we focus on two-partial 
incremental waves in this Section, we may take advantage of formulas established by 
Merodio and Ogden (2002) in a similar context. The components of the deformation 
gradient tensor F and of the left Cauchy-Green strain tensor B for 2D pre-strain and 
2D incremental motions, in the {x\,x 2 ) coordinate system (aligned with the (Xi,X 2 ) 
system), are 



K 
1 



B = 



1 + K 2 
K 



K 
1 



(4.6) 
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For plane strain, A 3 = 1, so that in incompressible solids A2 = X 1 1 by (2.1). It follows 
by (2.4) that I\ = I 2 . Accordingly, we define the single- variable function W{Ii) by the 
identity 

W{h) = W{h,h). (4.7) 
Then the 2D version of the constitutive equation (2.5) is 

er = -f,I + 2WiB, (4.8) 

where p is the Lagrange multiplier due to the incompressibility constraint. Also, the 
components of Aq are (Merodio and Ogden, 2002), 

Aojilk = ZWrfikBjt + AW lx B l3 B lk , (4.9) 

where T?i = W'(h), Wu = W"(h). With the help of (4.6) 2 , we find that in the (£i,£ 2 ) 
coordinate system, the non-zero components Aojuk relevant to in-plane motion are 

ioni2 - ioi2ii - 4WuK(l + K 2 ), 

Aim = A2111 - W\K + AW n K(l + K 2 ), 

A1122 = A)22ii = 4W n (l + K 2 ), .A01221 = A)2ii2 = AW n K 2 , 

-4oi2i2 = 2W\(1 + K 2 ) + AW n K 2 , A a2121 = 2W X + AW lx K 2 , 

-4 1222 = -402212 = 2WiK + 4W n K, .402122 = -402221 = 4W n K, 

iomi = 2W 1 (1 + K 2 )+AW 11 (1 + K 2 ) 2 , A Q2222 = 2t?i +AW U . (4.10) 

Now the remaining incremental governing equations (2.10)1,2,4 reduce to 

-Sii+iSi^-pw 2 ^, -S 12 + iS 22 = -pv 2 U 2 , iC/! + ^ = 0, (4.11) 
where by (2.7), 

iSn = i(ioim + + U01112C/2 + A 0U21 U[ + A 01122 U 2 - LP, 

LSi 2 = i4oni2f/i + i-4oi2i2?72 + (A1221 + p)^ + Aq\ 222 U 2 . 

\S 2 \ = L4oii2i^i + KA)i22i +p)U 2 + A 02121 U[ + A 02122 U 2 , 

\S 22 = i4oii22f/i + L4oi222?72 + ^02122^ + (A12222 + p)U 2 - LP, (4-12) 

where P is the increment of p. The pre-stress necessary to maintain the solid in the 
static state of large simple shear is (4.8). In particular, the <j 22 component along the x 2 
axis is found using (4.6) and (4.7) as 

a 22 = -p + 2W 1 , (4.13) 

leading to the connection 

^01221 +P = A2121 - 022- (4-14) 
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Note that Connor and Ogdcn (1995) and Fu (2005a) keep a 2 rather than 022 as a measure 
of the pre-stress. It is the component of the pre-stress along the principal axis x 2 , whose 
orientation changes with the pre-strain, in contrast to 022, the component of the Cauchy 
pre-stress tensor along the unchanged normal to the interface. As pointed out by Hussain 
and Ogden (2000), it is 622 which is continuous across the interface of two bonded sheared 
solids. 

Following the same procedure as in Section 3.1, we find the Stroh formulation of the 



governing equations in the expected form (2.10), where £ 
and N 2 , are given by 



[^1,^2, fti,^]* and -N, 



K 1 - 



022 



2Wi + 4WuK 2 




2W! 



AWuK 2 




(4.15) 



respectively, and — iV 3 by 

2(4W ? i- < 7 22 ) 



-(4W 1 -a 2 2)K 



-(AW! - a 22 )K 2W 1 K 2 + 2a 22 - 



^2 



2W 1 +4W 11 K 2 ] 



(4.16) 



4.2 Resolution 

The propagation condition (2.16) reduces to a quartic in q, 



2i 2 

pv ]q 



2(W 1 + 2W n K 2 )q A + 4(t?i + 2W n K 2 )Kq 3 

+ [2(Wi(2 + K 2 ) - AW n K 2 (2 - K 2 ) 
+ 4(Wi - 2WuK 2 )Kq 

+ 2V?i(l + K 2 ) + 4W n K 2 



pv 



0. (4.17) 



Notice how a 22 , though present in AT, does not appear explicitly in this equation. 

Quite surprisingly, there are two instances - both pointed out by Connor and Ogdcn 
(1995, 1996) - where we can solve this quartic exactly and simply. The first instance 
occurs for incremental deformations, when v — 0; it is applicable whatever the strain 
energy density W might be. The reason for the simplicity of this resolution is made 
apparent by the change of unknown from q to q = q — K/2. Then the quartic becomes 



2(Wi + 2W 11 K 2 )q i + [Wi(4 - K 2 ) - 2WnK 2 {4 + K 2 ) - pv 2 }? 

+ (pv 2 )Kq+ [{Wi + 2WnK 2 ){4 + K 2 ) - 2pv 2 ]{4 + K 2 )/8 = 0, (4.18) 

which is clearly a biquadratic at v = 0. The consequence is that any incremental static 
problem can be solved in its entirety for sheared solids because the roots q are accessible 
explicitly. This was to be expected though, because of an important theorem by Fu and 
Mielke (2002) which states that "the buckling condition for a pre-stressed elastic half- 
space is independent of the orientation of the surface as long as the surface normal remains 
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in the (£1,2:2) plane"; the buckling condition is what corresponds to the marginally 
stable static solution obtained at v = 0; it is also called the wrinkling condition, or the 
bifurcation criterion, or any other denomination associated with the onset of instability 
in the linearised (incremental) theory. 

The second instance where the quartic is easy to solve is when the solid is a Mooney- 
Rivlin material, see (3.22); this case is treated in Section 4.3. 

In general however, the quartic is difficult (but not impossible, see Fu (2005a, b), 
Destrade and Fu (2006), and the concluding Section) to solve analytically and other 
methods, such as those relying on the fundamental equations (2.25), are required. For 
the time being, we complete the picture with formal calculations. 

Assuming the roots of the quartic have been computed, and calling them qi, q 2 , q 1 , 
q 2 , where qi and q 2 both satisfy the decay condition (2.17), we find that the eigenvectors 
^ 1 and C 2 associated with q\ and q 2 , respectively, are 



(4.19) 



where 



& = [ij (721 {Oj + Kq 3 - 1) + (T22) , -(721 - £22)^ + J>i29j + 712 - pv 



21* 



(4.20) 



Here, 



721 = 2Wi + 4W U K 2 , 

$12 = Wx(2 + K 2 ) - 2W n K 2 (2 - K 2 ) 



7l2 = 2VFi(l + K 2 ) + AWnK 2 , 
z>i2 = 2WiK -AW n K 3 . 



(4.21) 



Then we compute the A = [a 1 ^ 2 ] and B = [b 1 ^ 2 ] matrices, and eventually the surface 
impedance matrix M 



as 



M = -i 



pv 



iBA 

721 (if + qi 
712 



52) 



72l(l + qiq 2 ) + (7 2 2 



- 721 - 022 {PV 2 - 712) 



qi + g2 
qiq2 



- Vvi 



(4.22) 



Note that this matrix is indeed Hermitian; this is easy to show by using the quartic (4.17) 
and the definitions (4.21) to uncover the identities: 



Qi + <?2 + Qi + q 2 = -2-K", 

qiq2q{q 2 = (712 - p« 2 )/72i, 

qm(qi + q 2 ) + qiq 2 (qi +92) = -2^12/721 ■ 

These identities allow us to rewrite the surface impedance matrix as 
72i((7i + qi - ~q~\ -q 2 )/2 72i(l + qiq 2 ) - a 22 

-721(1 + q{q 2 ) + CT22 -721 [qi 52(91 + q 2 ) -q{q 2 (qi +<72)]/2 



M = -i 



(4.23) 



(4.24) 
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Fu (2005a) derived the same form of the surface impedance matrix for the more general 
case of any plane strain (any pre-strain where A3 = 1), which includes the present case 
of finite shear. 

Notice also how the pre-stress CT22 is going to be explicitly present in the secular 
equation for Rayleigh surface waves (2.22), but absent from the secular equation for 
Stoneley interface waves (2.23). These secular equations remain implicit as long as the 
roots qi and q 2 are not known. In the general case where the quartic (4.17) is not solvable 
in a simple manner, we seek an explicit secular equation using the fundamental equations 
(2.25). 



Rayleigh surface waves. With the explicit expressions (4.15) and (4.16) for the blocks 
of the matrix N, we can compute TV -1 and N 2 . The lower left corner of these gives 
in turn K {1) = N 3 + pv 2 I and K {2 \ respectively. Then the equations (2.30) 

written at n = 1,-1,2, yield the linear homogeneous system 



(4.25) 



The vanishing of the determinant of the 3x3 matrix on the left hand side is the explicit 
polynomial secular equation for surface waves in a sheared incompressible semi-infinite 
solid. It is a polynomial of degree 4 in pv 2 . It is too long to reproduce in general but 
easy to obtain (and solve numerically) with a computer algebra system. Here we present 
its expression in the case where 022 = 0. Then, and K ^ are given by 



^11 


K (1) 
-^12 


K {1) 1 


K { - 1] 
^11 


^12 


-"22 


^11 


^12 


K {2) 
^22 _ 



E/i(0)Ci(0) 




"0" 


[/i(0)[/ 2 (0) + [/i(0)t/ 2 (0) 







U 2 {0)U 2 {0) 








pv 2 - 8Wx 4Wi 
AWi pv 2 + 2Wi 



2{4W 1 ~pv 2 )K 2W 1 {4-K 2 )-2pv 2 
2Wi (4 - K 2 ) - 2pv 2 -8W1K 



(4.26) 



respectively, the components of K ^ 1 ' are 



L 1 -2(2W 1 -pv 2 ), 



K 



(-1) 
12 



-2[2Wi(l + K 2 ) - pv 2 ]K, 
,,2 



= 4t? l( 2 + K 2 ) 2 + ^^-^(5 + ^-4^(1 + ^) 

W 1 +2W ll K 2 

(up to an inessential common factor), and the secular equation is the quartic 



(4.27) 



x 4 - 5x 3 + 




4 + K 2 



+ K 



2 Wu 



8 4- 



K 2 



K 2 



K 



■ Wit 
Wi 



W 1 + 2Wi 1 K 2 
Wx(A + K 2 ) 



0, (4.28) 



where x is the following non-dimensional measure of the squared wave speed, x = 
pv 2 /[W 1 (A + K 2 )]. 
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As stated above, the secular equation is also a quartic in the squared wave speed 
when i722 7^ 0. For a given material, a given pre-stress, and a given shear, its numerical 
resolution may yield more than one positive real root. If such is the case, then for 
each corresponding speed, compute the roots to the quartic (4.17) and discard those 
(supersonic) speeds which do not give two complex conjugate pairs of roots. Finally, find 
which of the remaining speeds (if there is more than one) satisfies the secular equation 
written in optimal form (2.22). 

Stoneley shear-twin interface waves. For Stoneley interface waves, the fundamen- 
tal equations (2.25) are not practical to derive secular equations in general, and we must 
resort to other methods, such as those developed by Destrade and Fu (2006). The ex- 
ception is the special case when each half-space is filled with Mooney-Rivlin materials, 
because then the roots to the quartic (4.17) can be found easily, see Section 4.3. 

We now focus on the possibility of propagating incremental waves at a shear-twin 
interface. In this configuration, two solids, made of the same incompressible material, are 
subject to equal and opposite shears, see Figure 5. The study of wave propagation at this 
type of interface can have important repercussions in the non-destructive evaluation of a 
twinned interface because this bimaterial can "simulate the finite (plastic) deformation 
associated with a crystal twin" (Hussain and Ogden, 2000) . 




Figure 5. Stoneley wave propagation at a shear-twin interface. Both half-spaces are 
occupied with the same incompressible solid, one subject to a simple shear of amount 
K, the other subject to a simple shear of amount —K. The analysis shows that Stoneley 
waves cannot actually travel in the direction of shear. 

For one half-space, the propagation condition is the quartic (4.17). For the other half- 
space, the amount of shear is changed to its opposite, but W\ and Wn do not change 
signs because they are functions of K 2 . It follows that if q% and q2 are qualifying roots 
of the propagation condition in one half-space, then —q\ and —q2 are qualifying roots in 
the other half-space. Also, we see from (4.21) that 721, 712, and $21 in one half-space 
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are equal to their counterparts in the other half-space, but that the counterpart to O12 is 
—ui2- Finally, &22 is continuous across the interface. Then we conclude from (4.20) that 
the counterparts to the matrices 



A21 A22 



B 



Bu B\2 
B 2 i B 2 2 



in one half-space are the matrices (Destrade, 2003; Ting, 2005) 



An 
-A21 



A12 

-A22 



B* = 



-B n 
B21 



— B\2 

B22 



(4.29) 



(4.30) 



in the other half-space. In turn, this conclusion leads to the following sum of surface 
impedance matrices 

~2M n 



M + M* = 







2M- 



22 



(4.31) 



where M* = iB*(A*) 1 . Consequently, the exact secular equation for Stoneley shear- 
twin interface waves is, according to (2.23), (4.24), and (4.31) that either 



3(<?i + q 2 ) = 0, or 3?[<7i<7 2 (<7i + q 2 ) 



0. 



(4.32) 



However, as is easily proved, neither of these quantities can be zero when both q\ and 
q2 have positive imaginary parts. It follows that Stoneley waves cannot propagate in 
the direction of sheared at a shear-twin interface, whatever the strain energy function is, 
and whatever what the shear is. Note however that Hussain and Ogden (2000) show, in 
their study of reflection and transmission of plane waves at shear- twin interface, that an 
incident harmonic plane wave can give rise to an interfacial wave (but of a different type 
than the Stoneley type). 

4.3 Examples 

Sheared Mooney-Rivlin solids. For the Mooney-Rivlin strain energy density (3.22) 
we have 



W 1 = V/2, W n = 0, 



(4.33) 



where V = T>\ +T> 2 is the shear modulus. The quartic propagation condition (4.17) then 
factorizes to 



(<7 2 + l)(q 2 + 2Kq + K 2 + rf ) = 0, where 77 = yjl - pv 2 /V , (4.34) 
and its roots with positive imaginary parts are 

qi = i, q 2 = ~K + wi. (4.35) 

We also have 

721=©, 7i2=^(l + ^ 2 ), £12 =V(2 + K 2 ), v 12 =VK, (4.36) 
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and the surface impedance matrix of (4.24) reduces to 



M = V 



i] + 1 
i(ry - 1) - \a 22 jV 



-K + i(rj-l)+i&2 2 /V 
rf +rj + K 2 



(4.37) 



For Rayleigh surface waves in a sheared Mooney-Rivlin material, the secular equation 
(2.22) is a cubic in rj, 



rj 3 + rj 2 + (3 + K 2 - 2a 22 /V)r, - (1 - a 22 /V) 2 



0. 



(4.38) 



See Connor and Ogden (1995) for this equation with a 2 instead of a 22l and Destrade and 
Ogden (2005) for a generalization of this equation to a triaxial stretch followed by a shear. 
See also Figure 6 for the variations of the squared scaled wave speed pv 2 /T> = 1 — rj 2 with 
the amount of shear K, for several values of the pre-stress a 22 . In particular, the plots 
at CT22 = ±21? show that the Mooney-Rivlin material is unstable when subject only to a 
hydrostatic pressure of that amount (that is, pv 2 ~ at K ~ 0) but regains stability as 
soon as it is sheared [pv 2 > at K ^ 0) . 



pv 2 /£> 

0.5 




Figure 6. Squared scaled surface wave speed pv 2 /T> as a function of the amount of shear 
K for a Mooney-Rivlin material. For the solid plot, there is no pre-stress normal to the 
boundary (a 22 = 0); for the dashed plots, the displayed number indicates the value of 
the non-dimensional pre-stress a 22 jT>. More plots are found in the article by Connor and 
Ogden (1995). 



As is clear from the roots (4.35), from the secular equation (4.38), and from Figure 6, 
both the penetration depth and the speed are even functions of K. Also, once rj is found 
by solving the secular equation, we have dr]/d(K 2 ) = —ri 2 /[2ri 3 + rj 2 + (1 — <t 22 /T>) 2 ] < 
and d^{q 2 ) / d(K 2 ) = drj/d(K 2 ) < 0, indicating that both the surface wave speed and 
the penetration depth increase with the magnitude of the shear. 
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We note that having found an explicit expression (4.37) for the surface impedance 
matrix allows us to solve the problem of Stoneley interface waves in its entirety, for half- 
spaces made of different Mooney-Rivlin solids, subject to different amounts of shear; see 
Chadwick and Jarvis (1979b) for more discussion on this point. 



Sheared Gent solids. As an example of a strain energy density for which the quartic 
(4.17) is not easily solved, we now work with the following function, 

W = -^CJ m hi(l-^^j , h<3 + J m , (4.39) 

proposed initially by Gent (1996) to describe strain-stiffening elastomers and since then 
extensively used by Horgan and Saccomandi (2003, and references therein) to model 
strain-stiffening soft biological tissues such as arteries. Here C > is the infinitesimal 
shear modulus and J m > is a constant, accounting for the limiting chain extensibility 
of the solid. Hence the condition (4.39)2 limits the amount of shear by imposing 

-^<k <yx; (4.40) 

Turning our attention to the propagation of surface waves, we compute the following 
quantities, 

_ C J c 7 

Wi = VFn = — - (4 41) 

To fix the ideas, we take two Gent materials, one with J m — 9.0, the other (stiffer) with 
J m = 1.0. Figure 7 shows the variations of the non-dimensional measure of the surface 
wave speed sj pv 2 /C with the amount of shear K . The bounds due to the inequalities 
(4.40) are clearly visible, indicating that the solids become more and more rigid as their 
limit of chain extensibility is approached. The speed is an even function of K and only 
the K ^ needs be displayed. For a boundary free of pre-stress (022 = 0), we solve 
numerically the quartic secular equation (4.28) and find in general that there is two real 
positive roots for pv 2 ; one gives a supersonic speed and is discarded; the other gives a 
speed for which the exact secular equation (2.22) is satisfied and it is kept. The Figure 
also shows (dotted plots) the effect of the compressive pre-stress 022 = —1.5C, which is to 
slow the wave down in the small-to-moderate shear region; here again a quartic secular 
equation is solved numerically and only one speed is kept. 



5 Propagation in a principal plane 

In this Section we consider an interface wave propagating in a principal plane, x-i = say, 
but not in a principal direction. We call 9 the angle between the propagation direction 
x\ and the principal axis x\, see Figure 8. Although some results exist for non-principal 
Stoneley waves (Destrade, 2005), the focus of this Section is on Rayleigh surface waves. 
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Figure 7. Scaled surface wave speed y/pv 2 /C as a function of the amount of shear K for 
two Gent materials. For the solid plots, there is no pre-stress normal to the boundary 
(<5"22 = 0); for the dashed plots, the pre-stress is 022 = — 1.5C. The vertical asymptotes 
correspond to the maximal amount of shear, that is K max = ±3.0 for the Gent solid with 



J m = 9.0, and K n 



±1.0 for the other, stiffer, Gent solid with J m = 1.0. 



5.1 Governing equations 

Hence we model this motion as 



{u,p, s} = {U(kx2), ikP(kx2), ikS(kx2)}e 



(5.1) 



where eg = cos 9, sg — sin8. By projecting the governing equations in the coordinate axis 
of the principal axes (where Aa has 15 independent non-zero components, see Section 
2.2), we can write them in the Stroh form (2.13) as a homogeneous linear system of six 
first-order differential equations, where £ = [U\, U2, U3, S21, S22, ^23]'; see Destrade et al. 
(2005) for details. Here the matrices —Ni, iY 2 , — AT 3 are given by 






ce{l2i - C2V721 





co 










Se(723 - <J 2 )/j23 






I/721 










X 





— K 
















// 











l/723_ 




— K 








(5.2) 
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Figure 8. Interface wave propagation when the boundary x 2 — is a principal plane of 
pre-strain. The x\ direction is the direction of propagation, making an angle 6 with the 
principal direction of pre-strain x\. 



respectively, where 

X = 2ce(/3i2 +721 - cr 2 ) + Sg73i, 

V = Cg[7i2 - (721 - Cr 2 ) 2 /7 2 i] + Sg[732 - (723 - CT 2 ) 2 /723], 

M = c e7i3 + 2s 2 , (#23 + 723-^2), 

k = c e s e {l3i3 - (3u - P23 ~ 721 - 723 + 2<7 2 ), (5.3) 

and the 7^, are defined in (3.12). 

Rogerson and Sandiford (1999) show that the propagation condition (2.16) is a cubic 
in q 2 , 

7217239 6 - [(721 + 723)^ - ci] q 4 + (X 2 - c 2 X + c 3 ) q 2 + (X - c 4 )(X - c 5 ) = 0, (5.4) 
with X = pv 2 and 

Cl = (721713 + 2/3l2723)cg + (723731 + 2/?2372l)s0, 

C2 = (723 + 7i3 + 2/3i 2 )cg + (721 + 731 + 2f3 2 3)sl, 
C3 = (712723 + 2/3i27i 3 )ce + (721732 + 2/?2373i)4 

+ [712721 + 713731 + 723732 - (/?13 ~ Pl2 ~ fe) 2 + 4/3i 2 /?23]c0Sg, 

C4 = 712 ce +732Se, 

C5 = 713C9 + 2/3i 3 Cgsl + 7ai4- ( 5 - 5 ) 
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5.2 Resolution for Rayleigh surface waves 

Finding the eigenvectors £ , £ 2 , £ 3 of N corresponding to the qualifying roots 91, q 2 , 
93 (say) is a lengthy task by hand, best left to a computer algebra system. In the end, 
we find that the C are of the forms: 



c 2 = 



(5.6) 



where 



ami + a 2Qi + ao 
-Qj + hq] + biq 3 
_d 4 qj + d 2 q 2 j + d _ 



Here the quantities m and n are given by 



hzq z 3 + hiqj 
[y - X)(q^ + mq? + n) 

gsqj + giq 3 



(5.7) 



rn = 



1 + %^4)[r,-X\ + (±+<»- 



721 721 { v - x ) 



— C2) 2 2 I 



723 



2k 



(721 ~ 0-2) (723 ~ cr 2 ) 
721723 - X) 



CgSg, (5.8 



x [( M -X)( ?? -X)-« 2 ]/( 72l723 ). (5.9) 

The expressions for the quantities at, hi, di, hi, gt are too lengthy to reproduce but they 
are easily obtained in a formal manner using a computer algebra system; as it turns out, 
these constants are not needed in the secular equation for Rayleigh surface waves. Indeed, 
the expression for the surface impedance matrix M is very lengthy, but its determinant 
factorizes greatly. Hence we find that the exact secular equation for Rayleigh surface 
waves (2.22) reduces to (Taziev, 1987; Destrade et al., 2005) 

wjj\ — u>in(m — (jJn) — 0, (5.10) 

where 

wi = -{qi + 92 + 93), wii = 9192 + 9293 + Q3Q1, w m = -919293- (5-11) 

This secular equation remains implicit as long as the roots qi, 92 , 93 satisfying the decay 
condition are not known. To compute them, we must find the wave speed, by using the 
fundamental equations (2.30). 

First, using the explicit expression (5.2) of the Stroh matrix N, we compute TV -1 
and AT 3 , and in particular we find explicit expressions for the lower left blocks K^ 1 ) = 
AT 3 + pv 2 I, K { ^\ K {3) . Next, we use the result (Destrade, 2005) that 17(0) is in the 
form 

U(0) = U 1 (0)[l,ia,(3] T , (5.12) 



" (721 ~ 02Y 
721 



4 + 



(723 - cr 2 y 



723 



■4 
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where a, (3 are real numbers, to write the fundamental equations (2.30) at n = — 1, 1,3 
as the following system of three equations, 



K-(-l) 

21 13 


^33 


^(-1)" 
-"■22 




"2/3" 








^(1) 
^33 






/? 2 






X (3) 


^33 


lx 22 




a 2 




__K-(3) 



Then, we solve this non-homogeneous system by Cramer's rule to find 

2/3 = Ai/A, /3 2 = A 2 /A, 
where A, A l7 and A 2 are the following determinants, 



(5.13) 



(5.14) 







^33 


A 22 


A = 


^13 


^33 


K {1) 

^22 




A 13 


A 33 


K (3) 
^22 



Ai 



^11 
^11 
^11 



K. 
K, 
K- 



(-i) 

33 
(1) 
33 
(3) 
33 



K. 



K. 
K. 



(-i) 

22 
(1) 
22 



(3) 



22 



A 2 = 



A 13 
^(1) 
^13 
K (3) 
^13 



-K { - 1] 
^Mi 

_if(3) 
^11 



Finally we write down the compatibility of equations (5.14) as 

4AA 2 = 0, 



A? 



K, 



(-i) 



22 

K {1) 

^22 
^22 



(5.15) 



(5.16) 



which is the explicit secular equation for non-principal surface waves in deformed incom- 
pressible materials. 

In general this equation is a polynomial of degree 12 in pv 2 (Tazicv, 1989; Destrade 
et al., 2005), easy to solve numerically. Of the 12 possible roots, we keep those which are 
real, positive, and give a subsonic speed (that is a speed for which the bicubic (5.4) has 
three pairs of complex conjugate roots). We then test the remaining speeds against the 
exact secular equation (5.10). 



5.3 Examples 

neo-Hookean solid. The neo-Hookean form of the strain energy density for an in- 
compressible isotropic solid is a sub-case of the Mooncy-Rivlin form, namely 2? 2 = in 
(3.22) so that 

W = 2?i(A 2 + A 2 + A 2 )/2. (5.17) 

It leads to a stress-strain relationship (2.5) which is "linear" with respect to the left 
Cauchy-green strain tensor, 

cr = -pI + V 1 B. (5.18) 
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The neo-Hookean model is unable to capture neither qualitatively nor quantitatively 
experimental data, over any range of deformations (Saccomandi, 2004). It is nonetheless a 
highly popular model in the literature because it forms the basis of a statistical treatment 
for the molecular description of rubber elasticity (Treloar, 1949). 

With respect to wave propagation, it has more peculiar properties than the Mooney- 
Rivlin material, because here, 



7«=©iA?, 2/3y = V 1 (X 2 + X 2 j ), 



(5.19) 



which leads to great simplifications in the matrices N\, N 2 , N% of (5.2). 

Now placing ourselves in the coordinates system {x\, £2, £3) attached to the wave 
propagation, sec Figure 8, we introduce the functions Ui, S 2 i (i = 1, 2, 3), defined by 



Ui = CtijUj, £>2i — ClijSzj, where flij = 



Cg -3g 

1 

Sg C 9 



(5.20) 



With these functions, the governing equations (2.13) decouple the anti-saggital motion 
[U3, S23] from its saggital counterpart (recall that the direction of propagation and the 
normal to the interface define what is called the saggital plane.) For the latter motion 
we find 

[Ui,U^J 21 ,S' 22 Y = iN[U 1 ,U 2 ,S 2 i,S22} t , (5.21) 



with 



N 





-1 
pv 2 - x 




-1 + C72 1/(©iA|) 0' 


0-1 



pv* — V 



-1 + CT2 







(5.22) 



where 



X = V 1 (c 2 X 2 1 +s 2 X 2 3 + 3X 2 2 -2W 2 ), u = V 1 [c 2 \\ + s 2 Xl + X 2 2 (l-n 2 ) 2 ] 1 (5.23) 

and (72 = o"2/(2?iA|) is a non-dimensional measure of the pre-stress. The associated 
propagation condition is 

(q 2 + l)[X 2 2 q 2 + (c 2 X\ + s 2 Xl) - pv 2 /V{\ = 0. (5.24) 

In other words, the situation is formally the same as that for principal waves in 
Mooney-Rivlin material, see Section 3.3. The conclusion is that the speed is given by 

pv 2 =V 1 (c 2 e X 2 + s 2 e Xl-X 2 2 r ) 2 ), (5.25) 

where 77 is the real root of (3.26). Flavin (1963) established this result at cr 2 = 0. As he 
also showed, the situation gets more complicated for Mooney-Rivlin solids, because the 
wave is no longer plane polarized for a triaxial pre-stretch. 
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Mooney-Rivlin solid. The Mooncy-Rivlin strain energy density (3.22) at T>2 ^ does 
not lead to a decoupling of the saggital motion from the anti-saggital motion. However, 
as in Sections 3 and 4, we find that the propagation condition factorizes, here as 



(q 2 + l)(q i -Sq 2 + P) = 0, (5.26) 

where 

\721 723/ \721 723/ \721 723/ 

P = (X - 7 12 ce - j 32 s 2 e)( x - 713C0 - 73iSe)/(72i723)- (5-27) 

Flavin (1963) was the first to notice that q\ = i is an attenuation factor for non- 
principal waves in Mooney-Rivlin solids. Pichugin (2001) showed that a necessary and 
sufficient condition for the factorization (5.26) to occur is that the relations (3.24) hold; 
he also showed, completing earlier work by Willson (1973a), that another factorization of 
the general bicubic (5.4) also occurs, this time for any strain energy function, when two of 
the principal stretches of pre-strain are equal (equi-biaxial pre-strain). Finally note that 
(q 2 + 1) always comes out as a factor in the propagation condition for inhomogeneous 
waves in Mooney-Rivlin solids, whatever the direction of propagation is (not necessarily 
in a principal plane as here), see Destrade (2002) for details. 

Thanks to the factorization (5.26), we can actually compute explicitly the quantities 
u>i, lju, Wni °f (5.11). Indeed, we now have 

qi = i, 9293 = -VP, q2 + q 3 = i^VP-S, (5.28) 

so that 

uti = -i(l + \J2VP-S), -oju = VP+ \J2VP-S, w„i = iVP, (5.29) 
leading to an explicit and exact form of the secular equation (5.10), 

n(l + \J 2VP - S 1 ) + Vp(m + VP + \J iVP - S 1 ) = 0. (5.30) 

This result allows us to investigate the influence of pre-stress on surface wave propaga- 
tion (see Destrade et al. (2005) for an example) and to address an important question in 
the study of surface stability: how much can a Mooney-Rivlin half-space be compressed 
before it buckles? In Section 3 we found the critical stretch in a principal direction (xi), 
indicating the appearance of wrinkles parallel to x 3 , see Figure 2 for a visualization. 
However, could it be that wrinkles appeared earlier in the compression, in another di- 
rection? To answer this we take X(= pv 2 ) = (onset of instability) and solve (5.30) for 
A, for each value of 9. 

For instance, take the case of the following plane pre-strain, 

Ai=A, A 2 =A-\ A 3 = l, (5.31) 
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imposed on a half-space made of the Mooney-Rivlin solid with material parameters 

Pi = 2.0/x, T> 2 = 0.8//, (5.32) 

where // has the dimension of a stiffness (the shear modulus of this Mooney-Rivlin solid is 
(T>i + T> 2 )/2 = 1.4//). Figure 9 displays the values of the critical stretch ratio, measured 
in the principal direction xi, for each angle 9 and for several values of the pre-stress a 2 - 
The solid line corresponds to a 2 — 0. At 9 = we find the critical compression ratio 
A" (say) for wrinkles aligned along x^; hence the solid curve starts at A° = 0.544 as 
expected from solving (3.27) with v — and Ai = A° = A^ 1 . At 9 ^ we find that the 
critical compression stretch is below A° ; it follows that A° is the absolute critical stretch 
of compression for our example (5.31), (5.32). Of course, for Mooney-Rivlin solids other 
than (5.32), or for pre-strains other than (5.31), or for solids other than Mooney-Rivlin 
solids, we might end up with a different behaviour in compression. However the same 
analysis can be brought in each case to its conclusion, with no additional difficulty. 



o.e - 
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Figure 9. Critical stretch ratio of compression for a Mooney-Rivlin material as a func- 
tion of the angle between the normal to the wrinkles and the principal axis of greatest 
compression. The Mooney-Rivlin solid is subject to a finite plane strain compression 
(A3 = 1); its material parameters are T>\ = 2.0/i, T>% = 0.8/j,; the pre-stress normal to the 
boundary is given by a 2 j [i = -3.0, -1.0, 0.0, 1.0, 3.0 (/i has the dimension of a stiffness). 



Gent solid. We find that the elastic moduli (3.12) for Gent materials (4.39) are given 
in general by 

Hi = c - 



2ft 



c 



A?- 



\ 2 

A 2 



X 2 



A,- . 



3 — Aj — A 2 



\ 2 



A? + \] + 



2(A? 



A?) 2 



Jn 



3 A-£ A2 A 



(5.33) 
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Now consider that a half-space made of a Gent material with J m — 9.0 is subject to 
a large shear such that the boundary is the plane of shear (in Section 4, the boundary 
was the glide plane) and the boundary xi = is free of tractions. Then the principal 
stretches are 

Ai = y/1 + K 2 /4 + K/2, A 2 = 1, A 3 = y/l + if 2 /4 - K/2, (5.34) 

and of course, \\ + \\ + A3 — 3 = K 2 . 

Figure 10 shows the variations of the surface wave speed in the plane of shear, with 
respect to the angle between the direction of propagation (x±) and the direction of shear 
(X±). For a shear of amount K = 0.5, the principal axis of greatest stretch x\ makes 
an angle 37.9° with the direction of shear, and the principal axis of smallest stretch X3 
makes an angle 127.9° with the direction of shear. For K = 1, those two angles are 31.7° 
and 121.7° , respectively; for K = 2, those two angles are 22.5° and 112.5° , respectively. 
Clearly, the surface wave reaches extremal values along those directions, indicating that 
the principal directions can be determined acoustically. 
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Figure 10. Wave propagating in the plane of shear of a semi-infinite sheared Gent solid 
with J m = 9.0: scaled speed pv 2 /C as a function of the angle between the propagation 
direction and the direction of shear, for several values of the amount of shear. 



6 Concluding remarks 

This Chapter has presented several situations where it is possible to derive explicit secular 
equations in exact form and in polynomial form, for incremental waves propagating 
along the plane interface of one (or two) deformed semi-infinite hyperelastic solid (rigidly 
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bonded solids). The subject of interface waves in general is broad and includes many 
other situations and geometries such as, to name but a few: the addition of a layer of finite 
thickness (Love waves, Lamb waves, etc.), or of a fluid (Scholte waves for inviscid fluids, 
Stoneley waves for viscous fluids, etc.), the consideration of an anisotropy due to families 
of extensible fibres, of cylindrical or spherical coordinates, or of curved boundaries, and 
so on. Reading the works by Ogden (2003, 2004) and by Guz (2002) and the references 
therein gives a good overview of the vast panorama spanned by these other types of 
interface wave problems. 

This concluding Section expands on the reasons put forth to explain why it can be 
advantageous at times to seek explicit secular equations rather than to turn to numerics 
outright. Among other things, explicit secular equations in polynomial form 
(i) are easy to solve numerically, with the greatest precision required; 
(ii) are sometimes surprisingly simple and short, in spite of a complicated or even 

unsolvable propagation condition; 
(in) lend themselves to simple asymptotic treatments, leading to approximate analytical 

expressions for the wave speed; 
(iv) account for all the solutions satisfying the propagation condition and the boundary 
condition at the interface. 

6.1 Numerics 

Point (i) goes without saying because the numerical methods used to determine the 
roots of a polynomial are safe and robust. Of course, at most only one root corresponds to 
the actual solution and all the other roots must be discarded as being spurious. So, once 
a likely root is found numerically (that is a root pv 2 which is real and positive), it must 
be checked that at that speed the exact secular equation is satisfied. This routine check 
is simple enough to perform (solve the propagation condition, find the corresponding 
impedance matrix, check that (2.22) (for Rayleigh waves) or that (2.23) (for Stoneley 
waves) is satisfied). 

Other numerical techniques to find the interface wave speed invoke the Barnctt-Lothe- 
Stroh integral formalism (see Ting (1996) for an account), computational algorithms for 
eigenvalue problems Taylor (1981), or an algebraic Ricatti equation for the impedance 
matrix (see the Chapter by Fu in this book and the references therein) to determine 
M{v) for any numerical value of v; then v is varied in order to satisfy the exact secular 
equation up to the desired precision. This however is not always an easy task, as seen in 
the following example. 

Consider the Mooney-Rivlin solid with material parameters (5.32), maintained in a 
state of static pre-strain, with (Rogerson and Sandiford, 1999), 

Ai = V31395, A 2 = V0J, A 3 = (AiA 2 ) _1 , a 2 = 0.8/i. (6.1) 

Take the interface between the solid and the vacuum to be the principal plane x 2 = 
and study the propagation of a Rayleigh surface wave in a direction close to X3 (9 is close 
to 90° in Figure 8.) 

Along X3 we have a principal wave, travelling with a speed v(90) say, found from 
Section 3.2. Here (Destrade et al., 2005) we find that v(90) ~ 1.327 y/fl/p. On the other 
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hand, a shear (homogeneous) bulk wave linearly polarized along x\ travels with speed v\ 
say, given by pv\ =731, and a shear (homogeneous) bulk wave linearly polarized along 
x 2 travels with speed v 2 say, given by pv\ — 732- Here we find that v\ ~ 1.384-\//z/p 
and that v 2 — 0.995^//i/p, showing that the surface wave travels with a speed which is 
intermediate between those of the shear bulk waves. That principal wave is two-partial 
and polarized in the (X2X3) plane, which is why it can afford to be faster than the shear 
bulk wave polarized in the x\ direction. Also, it is isolated, in the sense that the transition 
toward a surface wave propagating in a direction 9 ^ 90° is abrupt, because this latter 
wave is tri-partial and must therefore travel with a speed v(9) say, which is less than the 
speed of any homogeneous bulk shear wave, in particular less than \fc4Jp and less than 
yj c$/p, where C4 and C5 are defined in (5.5). Figure 11 shows the variations of these 3 
speeds in the neighbourhood of the X3 axis and makes those features clear. 




principal surface wave 

bulk shear waves 
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Figure 11. Wave propagation in the x 2 principal plane of semi-infinite Mooney-Rivlin 
solid (5.32) deformed by (6.1): variations of scaled speed y/pv 2 /p as a function of the 
angle, near the X3 direction. 

Now if our numerical method for the surface wave speed relies on using the known 
principal wave speed as an "initial guess" , then it might run into difficulties here because 
of the isolation of the principal wave. This special feature is characteristic of strongly 
anisotropic crystals in linear elasticity. In incremental dynamics, it can occur at will, 
simply by deforming the solid sufficiently to create a strong strain-induced anisotropy. 

Also, as is clear from the Figure, the surface wave travels with a speed which is 
extremely close to that of the bulk shear wave */ c^/p. Hence at 8 = 89.99°, the speed 
of the former is given by 

yjp~v r fp = 0.9948641879596778 ... , (6.2) 
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whilst the speed of the latter is given by 

= 0.9948641879596860 . 



(6.3) 



Consequently, if a numerical scheme for the surface wave speed relics on increasing v in 
small steps until det M(v) = is satisfied to any desired precision, then it will have to 
be pushed to the 14th significant digit to insure that the speed does not correspond to a 
bulk wave! Here the explicit polynomial secular gives only two positive roots, and it is a 
routine check to verify that only (6.2) satisfies the exact secular equation. 

6.2 Polynomials 

For an interface wave propagating in any direction in a deformed solid, the propaga- 
tion condition is a sextic in general. Such is for instance the case for a wave propagating 
in any direction in the glide plane of the sheared block in Figure 4. Although a sextic is 
unsolvable analytically, it is nonetheless possible to find a polynomial secular equation 
for a Raylcigh surface wave, see Taziev (1989) or Ting (2004) for details. The polyno- 
mial turns out to be of degree 27 in pv 2 . For Stoncley interface waves it does not seem 
practical to look for a polynomial secular equation. 

For an interface wave propagating in any direction in a principal plane of pre-defor- 
mation, the propagation condition is the bicubic (5.4). Although formulas exist for the 
roots of a cubic, the analytical resolution of (5.4) proves to be complicated because it 
is not known a priori whether the roots are purely imaginary or have a non-zero real 
part and accordingly, which formula should be used. However, the polynomial explicit 
secular equation for a Rayleigh surface wave is of degree 12 in pv 2 as seen earlier. This 
is also the case in the symmetry plane of an orthorhombic or monoclinic crystal. In 
the symmetry plane of a cubic crystal, the degree of the polynomial comes down to 10, 
as noted by Taylor (1981). For a two-partial surface wave, coupled to electrical fields 
through piezoelectricity in 2mm crystals, the polynomial is also of degree 10 in pv 2 for 
metallized boundary conditions (there the propagation condition is also a bicubic, see 
Collet and Destrade (2005)). 

For a two-partial (non-principal) interface wave polarized in a principal plane of prc- 
deformation, the propagation condition is a quartic, 

q 4 + d 3 q 3 + d 2 q 2 + d iq + d = 0, (6.4) 

say (see (4.17) for the case of simple shear). In contrast to the case of the bicubic above, 
we know here what the form of the roots should be, because the roots of a quartic are 
cither two pairs of complex conjugate numbers, or one double real root with one pair of 
complex conjugate numbers, or four real roots. Here only the first scenario is acceptable 
in order to satisfy the decay condition (2.17) and thus a single formula, always valid, is 
required; it can be found in the textbooks. First introduce in turn the quantities r, s, 
and h, defined by 

r = d 2 -fd§, s = dx - ±d 2 d 3 + \a%, h = d - \d x d 3 + j^d 2 dj - Jgrff, (6.5) 
the quantities A and <j> G]0,7r/3], defined by 

A= ^{I2h + r 2 f 2 , <p= i a rccos[^ (|.r 3 + s 2 - §r/i)] , (6.6) 
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and the quantities z\, z 2l and z 3 , defined by 



zi = 2A 1/3 cos(<?i>) - |r, 
z 2 = 2A 1/3 cos(0 + 27r/3) 
z 3 = 2A 1/3 cos(</- + 47r/3) 




(6.7) 



Then the qualifying roots are 



pi = isign(s)yzr- \d 3 + ±i(y/-z 2 + 
P2 = -Tjsign(s)^/IT - jd 3 + \\{^-z 2 



(6.8) 



where sign(s) equals 1 if s is non- negative and —1 otherwise. These formulas are per- 
fectly well handled by a computer algebra system, so that M, and ultimately the exact 
secular equation, can be found. It then means that any interface wave problem can be 
solved exactly. Nonetheless it might still be rewarding to look for the polynomial secular 
equation, to check whether they turn out to be simple. For instance, the polynomial 
secular equation for Rayleigh waves is a quartic in the squared speed when the solid is 
sheared (Section 4.2), or tri-axially stretched and then sheared (Destrade and Ogden, 
2005), or when the wave is polarized in the symmetry plane of a crystal (Currie, 1979). 
Now consider a flexural wave travelling along the free edge of a thin (Love-Kirchhoff) 
orthotropic plate where the edge makes an arbitrary angle with a principal axis of sym- 
metry. Thompson et al. (2002) show that the corresponding propagation condition is 
a quartic, but turn to a numerical resolution. However Fu (2003) shows that the poly- 
nomial secular equation is just a cubic in the squared speed. Finally consider the case 
of a piezoacoustic Bleustein-Gulyaev surface wave in a rotated F-cut about the Z axis 
for 4 crystals. There also the propagation condition is a quartic, but the fundamental 
equations reveal that the polynomial secular equation for metallized boundary conditions 
is simply a quadratic in the squared speed (Collet and Destrade, 2004). Clearly it is a 
worthy enterprise to unearth these polynomials rather than use numerics or the formulas 



6.3 Approximate expressions 

Once the polynomial secular equation is established, in the form V(v) = 0, say, it is a 
straightforward matter to derive an analytical approximation for v. Calling v n an initial 
approximation, we find in the first order that 



Of course we must choose vq judiciously for an optimal expression. 

Here we work out a simple example. We take a solid half-space subject to a hydrostatic 
pressure only, so that Ai = A 2 = A 3 = A say, and o\ = a 2 = a 3 = a say. From 
the incompressibility constraint (2.1), A = 1 follows and we have a pre-stressed, but 
unstrained, solid. It is thus isotropic and a surface wave propagates with the same speed 
in every direction. To derive the secular equation, we specialize for instance (3.20) to the 



(6.5)-(6.8). 



v ~ v - V(vo)/V'(vo). 



(6.9) 
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isotropic case. We find that 712 = 721 = /?i2 = Mch the infinitesimal shear modulus, and 
the exact secular equation is (Dowaikh and Ogden, 1990) 



f( V ) = V 3 + V 2 + (3 - 2a) V - (1 - of = 0, 



(6.10) 



where a = a / /irj and rj = y/l — pv 2 / hq. Calling c = y/ pv 2 / /irj a dimensionless measure of 
the wave speed, and multiplying (6.10) by /(—??), gives the polynomial secular equation 
(Dowaikh and Ogden, 1990), 



T(v) 



4(2 - a)c 4 + 6(ff - 2) 2 c 2 + (a 2 - 4)(a - 2) 2 



0. 



In the region where c is close to 1, the first order approximation (6.9) gives c 
V(l)/V'(l), that is 



21 - 28a + 6a 2 + 4ct 3 - a 4 



(6.11) 



(6.12) 



22 - 32a + 12a 2 

In particular we find that for an unstressed incompressible solid, c ~ 21/22 = 0.9545 . . ., 
which is a good approximation to the value found from the exact equation 0.9553. . . (Lord 
Rayleigh, 1885). On Figure 12 we superpose the approximate and exact curves, and find 
good agreement in the cr = — 1 to (7 = 1 region. Beyond that range, c departs from the 
neighbourhood of 1, and a better choice must be made for the initial guess. This point 
is not developed further here but it is clearly exposed in the paper by Mozhaev (1991). 
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Figure 12. Surface wave speed in an incompressible solid subject to hydrostatic pressure only. 
Solid plot: exact calculation, dashed plot: approximate calculation. 



6.4 Other waves 

Taylor (1981) and Taziev (1989) both remarked that many of the roots to an explicit 
secular equation in polynomial form have a physical origin, because they correspond 
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to a motion satisfying both the equations of motion and the boundary conditions at 
the interface. In particular for the solid/vacuum interface, the polynomial provides the 
parameters not only for the Raylcigh surface wave, but also for pseudo-surface waves, 
for Brewster reflection of bulk waves, for the reflection of inhomogeneous waves, for 
"organ-pipe" modes, etc. 

One of the most important root is that corresponding to a leaky surface wave (if it 
exists), whose energy is diffused slowly into the half-space, and for which the classical 
methods evoked in Section 6.1 run into major difficulties. 

Another application for this wealth of information is that many more wavefronts can 
be found, in particular those observed in the neighbourhood of cusps, which correspond 
to the interference of inhomogeneous plane waves (Huet, 2006). Figure 13 shows the 
remarkable correspondence obtained between experiment and predictions, for wavefronts 
propagating on the surface of a copper sample (cubic symmetry) . Here the explicit secular 
equation proves useful because it is easy to differentiate, as required for the derivation 
of the curve. 




Figure 13. Wavefronts on the surface of copper cut along a symmetry plane. Left: experi- 
mental results obtained by LASER impact (Huet, 2006). Right: predictions derived from the 
secular equation in polynomial form, showing in black, the Rayleigh wave wavefront and in grey, 
the wavefronts due to leaky waves and interference of inhomogeneous waves. 
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